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Effects of Geometrical Symmetry on the Vortex Nucleation and Penetration in 
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We investigate how the geometrical symmetry affects the penetration and arrangement of vor- 
tices in mesoscopic superconductors using self-consistent Bogoliubov-de Gennes equations. We find 
that the entrance of the vortex happens when the current density at the hot spots reaches the de- 
pairing current density. Through determining the spatial distribution of hot spots, the geometrical 
symmetry of the superconducting sample influences the nucleation and entrance of vortices. Our 
results propose one possible experimental approach to control and manipulate the quantum states 
of mesoscopic superconductors with their topological geometries, and they can be easily generalized 
to the confined superfiuids and Bose- Einstein condensates. 
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The advances in modern nanotechnology and the de- 
velopment of quantum computing have opened many new 
perspectives for research on mesoscopic superconductors 
One fascinating aspect in mesoscopic superconduc- 
tivity is the novel physics associated with vortices, vifhich 
has been a subject of great experimental [3-[loj and the- 
oretical interest in the past decades. In the meso- 
scopic superconductors, vortices are quantized and con- 
fined by the sample geometry, so they can exhibit many 
unique phenomena compared with conventional bulk su- 
perconductors. For example, when the sample is mesosp- 
copic, the giant vortex can form, and there may exhibit 
the exotic paramagnetic Meissner effect [l^. Further, 
it is shown that the symmetry of the sample geome- 
try can dramatically affect the properties of the meso- 
scopic system Recent experiment has also reported 
the symmetry- induced antivorticcs formation 0, which 
is a unique character of a mesoscopic superconductor. 
However, how the geometrical symmetry of sample can 
control the vortex state is hitherto not well understood. 

In this Letter, we explore the effect of symmetry on 
the vortex nucleation and entrance in the mesoscoopic 
superconductor far below the critical temperature. We 
develop an effective numerical method to generally solve 
the Bogoliubov-de Gennes (BdG) equations 1151 . w hich 
is based on the finite element method (FEM) [lil . [l7j . 
Compared with the conventional Ginzburg-Landau (GL) 
theory [18|, the BdG equations work well in wider tem- 



perature region, and can give the spectrum of excita- 
tions for spatially inhomogeneous superconductor sclf- 
consistently, supposed the few fundamental material pa- 
rameters are given. Several groups have successfully used 
the BdG equations to study the single vortex line [19|-|22[ . 

Here, we solve the BdG equations self-consistently, 
and obtain the vortex pattern and the current density 
in mesoscopic superconductors with arbitrary and com- 



plicated geometries. Our results show that the current 
distribution and the penetration of vortices are deter- 
mined by the symmetry of the sample geometry, and the 
entrance of the vortex happens only when the current 
density at the hot spots (i.e. the spots with maximum 
current density) reaches the depairing current density. 
These facts reveal unambiguously that the geometrical 
symmetry of the superconducting sample influences the 
nucleation and entrance of vortices, through determining 
the spatial distribution of hot spots. These results provide 
a practicable route to manipulate the quantum states of 
a mesoscopic superconductors in future applications. 

We start with the BdG equations for the quasiparti- 
cle wave functions u„(r) and Vn(r) in the presence of a 
magnetic field 
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Un{v) + A(r)w„(r) = £;„u„(r), 

w„(r) + A*(r)u„(r) = £;„u„(r), 

(1) 



where £'„ is the n-th energy eigenvalue, A(r) the pair 
potential, A(r) the vector potential and ^ the chemical 
potential. 

The pair potential A(r) is determined self-consistently 

by 

n 

A(r)=.g "n(rK(r)(l-2/(i?„)), (2) 

\En\<Ea 

where g is the interaction constant, f{E) the Fermi distri- 
bution function and Ec the cutoff energy which is related 
by the BCS relation via the transition temperature Tc 
and the superconducting energy gap Aq. 
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The vector potential A(r) is related to the current dis- 
tribution j(r) by Maxwell's equation 



Air 

V X V X A(r) = — j(r), 
c 



(3) 



where the current distribution [151 . |21| is given by 



/(ii;„K(r)(V--A(r)K(r) 



+ (1 - fiE,Mr)iV -A(r))<(r) - H.c. 
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The chemical potential /i is determined by the particle 
number conservation imposed on this system 

^ = 2/E {f{En)\u^{v)? + (1 - 2!{EMvn{v)?) ,(5) 



where N is the total number of particles. 

The boundary conditions for the above equations are 
given by 



h e \ / ft e , 

n • I -V - -A u„ = 0, n • -V + -A w„ = 0, (6) 



where we consider two dimensional superconducting sam- 
ples placed in the (x, y) plane, which are immersed in in- 
sulating medium in the presence of a perpendicular uni- 
form magnetic field along z direction, n is the normal 
vector of the boundary. 

We numerically solve Eqs. (|T|-(|6t s elf-consistently, 
based on finite elements method [ly, |l7[ • AH input pa- 
rameters used in the calculation are microscopic parame- 
ters which can in prin ciple be obtained from band struc- 
ture calculations 121[. These parameters consist of the 
cutoff energy E^, the coupling constant g and kp^Q [loj . 
where kp (vp) is the Fermi wave number (velocity) and 
^0 ~ fivp/Ag the coherence length. 

In this paper, we consider the system at the temper- 
ature T = O.lTc and choose E^ = 5Ao and kp^Q = 2. 
To conveniently compare our results with the experi- 
ment, other parameters are chosen so as to make the 
Ginzburg-Landau parameter k to be a specified value, 
where k = 0.96Al(0)/5o and Ai,(0) is the London 
penetration depth. We consider a square superconductor 
with the length of a side a = 5^o and k = 20. The sam- 
ple is first cooled down into the superconducting state, 
and then we slowly increase the magnetic field and solve 
the BdG equations at each field, mimicing the zero field 
cooling measurements [l, [3| • The global magnetization 
M can be calculated as M = 1/(25) r x jdxdy, where 
S is the sample area. In the frame of BdG theory, the 
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FIG. 1: (Color online) Magnetization (M) versus the external 
magnetic field (H) for a superconducting square sample with- 
out (thin line) and with (thick line) one defect. The vorticity 
L are labeled near the branches. For the sample with one 
defect, the corresponding vortex pattern at each L is given as 
the contour plots on the insets. 



free energy T of the system is given by 
^ = 2 ^ EnfiEn) -2j2En f rf'r|w„(r)|2 
^2kBTY,{.f{En)lnf{En) + [1 - .f(En)]\n[l - f{E,,)]} 



,3,^^/,3,|H(r)-HoP^ 

Stt 



where kg is the Boltzmann constant, Hq the external 
magnetic field and H(r) the spatial dependent magnetic 
field inside the sample. 

Shown in Fig. 1 is the Hq dependent magnetization 
of the square sample. We can see that the magnetiza- 
tion shows a zigzag-like pattern. Each descending branch 
of the zigzag can be described by a fiuxoid number L 
(vorticity), which determines how many times the phase 
of pair potential A(r) change by 27r along the sample's 
circumference. Here L is just the number of vortices. 
Sweeping the magnetic field continuously, the magneti- 
zation evolves along one of the fiuxoid curves until it 
reaches its end and jumps to the next curve, belonging 
to another fiuxoid state. 

In a perfect square sample with four-fold symmetry, 
we find that vortices enter the sample in unit of 4 with 
L = An (n = 0,1,2,...), and they replicate the geomet- 
rical symmetry of the sample, see Fig. 1. Note that the 
vortices enter the sample as 4 individual single quantum 
vortices rather than a giant vortex, which is in agreement 
with the previous works |17l. l24l|. However, if we reduce 
the sample size or increase the temperature, the vortices 
merge into one giant vortex. Further, upon the entrance 
of vortices, the magnetization displays a sudden jump to 
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the next branch. And correspondingly the free energy 
of the system drops discontinuously. This is consistent 
with the fact that the superconductor under considera- 
tion here is a type II one, which characterizes a negative 
surface energy Besides, we also study samples with 
other geometries, such as rectangle which has a two-fold 
symmetry, and disk which is rotation-invariant. In the 
case of rectangle, our results show that the vortices pen- 
etrate the sample in pair, just as expected. While for 
a disk, the Meissner state persists until H > 12<I>o/S', 
indicating that the critical field for vortex penetration 
is much larger than other geometry, which has been re- 
ported in experiment Q. Thus, the above results suggest 
that the geometrical symmetry of sample determines how 
the vortices entre the sample. 

When the symmetry of sample is broken by placing one 
defect on the sample boundary, the situation is dramat- 
ically changed. As an example, we consider the square 
sample with one sharp defect located at 2/5 of the length 
at the bottom edge, which mimics the practical experi- 
ment realized in ref. While in the perfect sample, the 
vortices enter in unit of 4, in the sample with one defect, 
the vortices enters one by one. Furthermore, from the 
vortex pattern (see the inset in Fig. 1), we can see that 
when the number of vortices inside the sample becomes 
large, the increased interaction between vortices domi- 
nates over the influence of the boundary, which makes 
the triangular lattice a favorable arrangement. The sys- 
tem free energy follows the same trend as in the per- 
fect square, once a vortex enters the sample, the free 
energy of system makes a discontinuous drop to a lower 
energy level. This confirms the determining influence of 
the sample symmetry on the vortex penetration and ar- 
rangement. 

Now we address how the geometrical symmetry deter- 
mines the vortex entrance and which intrinsic material 
parameter controls the nucleation and entrance of the 
vortices. We calculated the current density inside the 
samples j(r). The current distributions at several rep- 
resentative fields are shown in Fig. 2. It is evident that 
the current distribution conforms to the geometrical sym- 
metry. In the perfect square sample (Fig. 2, contour 
plot I), the current density distribution has a four-fold 
symmetry, forming 4 equivalent hot spots where the cur- 
rent density is the highest. Upon increasing the mag- 
netic field, the current density at the hot spots increases 
monotonously until it reaches a critical current jc simul- 
taneously. Further increasing the magnetic field results 
in the nucleation of vortices at the hot spots, and at 
the same time the current density drops to a magnitude 
much smaller than jc- For the sample with one defect 
(Fig. 2, contour plot II), the process is the same except 
that there is only one hot spot, due to the addition of the 
defect. Consequently, the change of current distribution 
leads to the large deviation of vortex entrance behavior. 
We notice that the critical current density for the vor- 
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FIG. 2: (Color online) The spatial dependent current density 
of the sample. The insets (I) and (II) are the contour plots 
of the magnitude of current density for the sample without 
and with a defect, respectively. The lines corresponds to the 
intersecting line of the current density in (I) and (II). The 
solid line a and d are for the perfect sample under the mag- 
netic field H = S-S^o/S (a, before the vortex entrance) and 
H = 8.4$o/'S' (d, after the vortex entrance); the thin fine b 
and c are the intersecting line across the defect, under the 
magnetic field H — 5.6$o/<S' (c, before the vortex entrance) 
and H = 5.8^o/S (b, after the vortex entrance); 

tex nucleation jc 2.5enh/m£^o is roughly equal to the 
depairing current density derived from London theory, 
Jc = cHc/iwX = 2.56enh/ni£^o. These facts suggest that 
in the mesoscopic superconductors, the intrinsic material 
parameter that determines the nucleation of vortices is 
the depairing current density. 

We further consider the cases where there are two de- 
fects on the sample boundary. When the defects dis- 
tribute symmetrically, for example, the two defects are 
at the middle positions of opposite edges, the sample 
has two fold symmetry and two equivalent hot spots are 
formed around the two defects. Thus the vortices enter 
the sample in pair through the defects, as shown in Fig. 
3a. On the contrary, if the defects are not equivalent, 
there is a difference between the maximum current den- 
sity at the defects. As a result, the vortex penetrates the 
sample one by one, and the vortex nucleates first at the 
defect with higher maximum current density. 

Therefore, our results support such a scenario for the 
vortex nucleation and entrance in the mesoscopic super- 
conductor: When an external magnetic field is applied, 
before the penetration of vortices, the magnetic field is 
screened out by the circulating current, resulting in a 
Meissner state. The spatial distribution of current is de- 
termined by the symmetry of the sample geometry and 
the distribution of defects. The points where the magni- 
tude of current are maximum will form hot spots. As the 
magnetic field increases, the current density at the hot 
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FIG. 3: (Color online) The contour plots of the magnitude of 
current density of the sample with two defects at the magnetic 
field H = 6.2$o/S. At this magnetic field, two vortices (left) 
and one vortex (right) has penetrated into the sample. The 
left is the sample with two defects arranged symmetrically at 
the middle of up and bottom edges, while in the right, two 
defects are not equivalent. 



spots increases correspondingly. Once the current den- 
sity at the hot spots reach the depairing current density, 
vortices will nucleate at the hot spots and enter the su- 
perconductor. Accompanying the entrance of vortex, the 
system magnetization and free energy have one discon- 
tinue drop and the system jumps to the next quantum 
state through a first order phase transition. Thus the 
geometric symmetry of the sample influences the spatial 
distribution of the hot spots, resulting in different vortex 
entrance behaviors. The influence of geometrical symme- 
try is a general property of the confined vortex system, 
thus this scenario can be easily generalized to the con- 
fined superfiuids and Bose-Einstein condensates. 

In conclusion, by developing an effective numerical 
method to solve the BdG equations self-consistently, the 
effect of the geometrical symmetry on vortex penetra- 
tion in mesoscopic superconductors and its mechanism 
are studied quantitatively. We demonstrated that, the 
condition of the nucleation of vortex is that the current 
density at hot spots reaches the depairing current den- 
sity. The geometrical symmetry influences the vortex 
nucleation and entrance through determining the spatial 
distribution of hot spots. The entrance of vortices leads 
to one first-order transition between the quantum states 
with different number of vortices. Our results suggest 
that by modifying the topological factors such as geomet- 
rical symmetry and defects etc., it is possible to control 
the entrance and location of the vortices, and thus the 
quantum states. In practical experiment, one can control 
the quantum states of a mesoscopic superconductor by 
tuning the applied magnetic field or current density. This 
opens up an applicable way to manipulate the quantum 
states in mesoscopic superconductors, which is crucial for 
nano-devices based on mesoscopic superconductors. 
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